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Abstract 

\C . 

, We study in a systematic way all static solutions of the Goldstone model in 

o 

1+1 dimension with a periodicity condition imposed on the spatial coordinate. 

OC . 

Q\ , The solutions are presented in terms of the standard trigonometric functions 

and of Jacobi elliptic functions. Their stability analysis is carried out, and the 
Oh. complete list of classically stable quasi-topological solitons is given. 
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1 Introduction 



The Goldstone model, the 0(2) invariant Higgs model with Mexican hat potential, 
is well known as a theoretical laboratory for the study of spontaneous symmetry 
breaking in relativistic field theory Recently, it was pointed out that the structure 
of the classical solutions of its 1+1 dimensional version on a spatial circle is strongly 
reminiscent of that of the localized solutions of the two-Higgs-doublet extension of 
the standard model (2HSM) of electroweak interactions |l|. Several branches of its 
solutions and their analogs in the 2HSM were explicitly constructed. Furthermore, 
the Goldstone model on S 1 is the simplest paradigm of field theories [|IJ, 0], [|J whose 
topological properties are too trivial to lead to absolutely stable solitons of any kind, 
and which nevertheless support the existence of classically stable quasi-topological 
solutions for some range of their parameters, here the Higgs masses or equivalently the 
radius L of spatial S 1 . The purpose of this little note is to present the complete list of 
static classical solutions of the Goldstone model on S l , the corresponding bifurcation 
tree and their stability properties. This analysis, apart from its own mathematical 
interest, may provide insight useful in our search for stable solitons in the 2HSM or 
in the special case of the minimal supersymmetric standard model. 

2 The classical solutions 

In this section we present all static classical solutions of the Goldstone model on 
spatial S 1 . The model is defined with two real Higgs fields (f)% and 02, whose dynamics 
is described by the action 



C 



^(0 1 ,0 2 ) = i(0? + ^-l) 2 



(1) 



(2) 
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and with periodic boundary condition on the space coordinate x G [0, 2irL]. The 
energy functional for static configurations is given by 



E = dx 



As usual, the static solutions of this model may be thought of as classical periodic 
motions of period 27rL of a particle in two dimensions under the influence of the 
inverted potential — V(0i, <p2)- Apart from the vacuum solutions, which have E vac = 
0, there is a trivial solution 

Ln 

1 = 2 = O , £ = — (4) 

which exists for all values of L. The corresponding small oscillation eigenmodes, 
labelled by j, have 

Co 2 (j) = ^(f-L 2 ) , j = 0,1,2,... . (5) 

They are four times degenerate, apart from the j = one which is doubly degenerate. 
This solution is always unstable since cD 2 (0) = —1/2. 

Many additional solutions, some of which were discussed in |I[ , bifurcate from 
the solution <fi\ = <p2 = at critical values of L. The generic solution of the model 
(|l|) can be expressed in terms of Jacobi elliptic functions; it has the form 

<pi = y/RcosQ , 02 = v^RsinO (6) 



R(x) = ai + a 2 sn 2 ( V~2A(x — x ), k 



nw =c L'w) dv (7) 

where a\, A, £, k, C, xq are constants, while sn denotes the Jacobi elliptic function; 
sn(z,k), sn(z,k) 2 are periodic functions on the real axis with periods 4K(k) and 
2K(k), respectively. K(k) is the complete elliptic integral of the first kind. Inserting 
(Q) and (0) into the equations corresponding to (0) leads to several conditions on the 
parameters and correspondingly to the following three types of non-trivial periodic 
solutions: 
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2.1 Type-I solution 

For a 2 = the function R(x) is a constant and the solution reduces to the form 



/ N 2 [Nx \ 1 N 2 . [Nx . 

1-— cos — + 6) , 2 = Wl- sm — + ^ (8) 



L 2 \ L J V £ 2 U 

where A/" is an integer and 9 an arbitrary constant; it was first obtained in ||. This 
solution bifurcates from ([|) at L = N i.e. when one of the uj of Eq. (|5|) crosses zero. 
The Higgs field winds N times around the top of the Mexican hat. Its energy is given 
by 

ttN 2 N 2 

E * L > ir > = -2L<?-J?) ( 9 ) 
2.2 Type-II solution 

For C = the function Q(x) vanishes, the parameter % vanishes as well and the 
solution takes the form 

1 = 2fcAsn(v / 2Ax,A;) , 2 = , A 2 = ^ — (10) 

The Higgs field oscillates in the <p2 = plane about the origin 0i = = <p2- Evidently, 
this solution is the Manton-Samols sphaleron [|J embedded into the Goldstone model. 
The argument k of the Jacobi elliptic function sn has to be chosen such that the 
periodicity condition is fullfilled, i.e. 

2K(k)m 
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-VTTk^ , (li) 



for some integer m. When k — > (i.e. L — > m) the solution flTOf) approaches (Q). 
Note that x can be translated by a constant and that the field + i<p2 can be rotated 
by a constant phase. The energy En is identical to the one of the one-Higgs model 
0. The relevant integral reads 



1m fK(k) 



l + k 2 ^ 2k 2 -l-k 4 - 



E n (L,m) = —!= / dy (k 2 sn 2 (y,k) ) + 12 

A ; a/2A(1 + k 2 ) 2 Jo LV K ' 2 > 8 J v ' 
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and can be evaluated in terms of the elliptic functions K(k), E(k) by means of the 
integrals 

K-E 



J K{k) dy sn 2 (y,k) 



k 2 



In particular 



TY17T 

for k = Eu(L = m, m) = — 

for k = l E n (L = oo,m) = -^t (14) 

3v 2 

2.3 Type-Ill solution 

When all the parameters entering in (Bp, (I7p are non zero the solutions are more 
involved. They were mentioned without details in || and we construct them here, as 
explicitely as possible. The equations imply the conditions 

ai = -(1-2A 2 (1 + A; 2 )) , a 2 = 4k 2 A 2 (15) 
3 

C 2 = — (l + (4k 2 - 2) A 2 ) (l + (4 - 2A; 2 ) A 2 ) (l - (2A; 2 + 2) A 2 ) (16) 

leaving a family of solutions depending on the four parameters : k, A, £, x$. The 
condition that R and C 2 should be positive implies 

A 2 < - 1 — (17) 
~ 2(1 + k 2 ) V ; 

When the equality holds one is led to the type-II solutions discussed above. 

In order for Q and R to be periodic on [0, 2irL] the following conditions must be 
satisfied 

C / -j^-dy = 2nn (18) 
Jo R(y) 

L = ^l (19) 

for some positive integers m, n. These conditions fix A and k as functions of L, m, n. 
Thus, for a given L, the generic solution depends on the parameters £, xq, m and 



n. The first two correspond to the arbitrary global phase and position of the con- 
figuration. The integer parameters m and n determine respectively the number of 
oscillations of the modulus of the Higgs field and the number of the Higgs field wind- 
ings around the origin <pi = = <p2 in a period 2ixL. 

Solving (18) , Ql9D in the case k = (with K(0) = vr/2) we find easily the critical 
values of L where the type-Ill solutions start to exist : 

A 2 = , "f => L 2 = \bn 2 -m 2 ) (20) 

4(bn 2 -m 2 ) 2 V ; v ; 

The expression for A 2 above combined with (17) lead to the condition 2n > m on the 
integers m and n. The positivity of L 2 in (20) then follows automatically. This result 
also demonstrates that for n fixed there are 2n — 1 possible values of m. 
Due to the absence of a closed form for the integral 

rK(k) i 

;dy (21) 



c + sn 2 (y, k) 

for generic values of k, the condition (|18| ) is impossible to handle algebraically. Even 
so, one can make some progress with the analysis by means of a k 2 expansion. 

For any n and m ^ 2n the coefficient of sn 2 in the integral giving Q(x) in (0) is 
proportional to k 2 and one may easily expand the solution in powers of k 2 . We find 

L 2 = e -?l^ + 0(k') 

^ 9 2n 2 (m 2 -4n 2 ) 2 . , . 

Inspection of the limit k = shows that the 2n — 1 solutions of type-Ill bifurcates 
from the type-I solution with N = n at L 2 = (6n 2 — m 2 )/2, for m = 1, 2, . . . 2n — 1. 
The energy can also be expanded in powers of k 2 leading to 

^ /, x 7i\^2n 2 / 9 9 . 3m 2 (12n 4 — m 2 n 2 — m 4 ) 

(ott 2 — m^)^ v bA{bn z — m z ) {An z — m z ) ' 
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The dependence on L is recovered by fl2"2"|). Correspondingly, the £; 2 -expansion of the 
energy of the Type-I solution about the point L 2 = (6n 2 — m 2 )/2 leads to 

ny/2n 2 (r .^ 2 _ 2 ;4 3m 2 (4n 2 + m 2 )(3n 2 - m 2 ) t /on6 , 



Ej(k, N = n)= (5n 2 - m 2 - -i + 

(24) 

The case m = n = 1 corresponds to the branch labelled W\ in []TJ . In this case the 
energies Ej and Em deviate only from the k 8 term on : 

1) = 4J[1 - ^ + f ) - + (*»)] (25) 
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E in (k, 1, 1) = 7T— \ - 1 (F + k b ) k* + C(A; iU ) (26) 

v ' ' ; 25V 2 L 128 V ; 294912 v ; J K ' 

Thus Ej(k, 1) is lower than Em(k, 1, 1) but only very slightly, as pointed out in [|TJ 
on the basis of a numerical study of these solutions. 

We have studied Eqs. (|i~8|), (O) numerically, solving for A 2 as a function of k 2 for 
different values of n/m. For fixed values of k, n/m we find a single solution A 2 (/c, n/m) 
obeying the following property 

7T) 2 1 

A 2 (fc = 0,n/m) = 4(6 ^ 2 _ m2) , A a (* = l,n/m) = - (27) 

This is illustrated in Fig. [1] for n/m = 1 and n/m = 4/7 by the solid lines, the dashed 
line representing the limit (|17|) . The second limit ( P7|) and the form of the solutions 
suggest that in the limit k — > 1, which corresponds to L — > oo, solution III aproaches 
solution II. To test this statement, it is interesting to compare their energies. The 
energy of the type-Ill solution is given by the integral 

rff^ pK(k) r , \ 2 1 / \n 

E m (k,m,n) = -^J Q ^[(ai-l+a 2 sn 2 (|/,A;)) +-(l + 16A 4 (A; 2 -l-A; 4 ))] (28) 
For k = 1 one obtains 

Emih m, n) = — ^ = -E n (l, m) (29) 

This verifies our expectation that solution III approaches solution II in the limit 
k — *■ 1. The occurence of the factor 1/2 is due to the fact that the solution III 



depends on sn 2 which has a period 2K(k) while the solution II depends on sn whose 
period is 4K(k). 

The energies of the solutions of a few low lying branches are plotted in Fig. [2] as 
functions of L. For L > y~5/2 all four types of solutions coexist and satisfy 

Ej(L, N = 1) < E H i(L, m = 1, n = 1) < E n (L, m = 1) < E (L) (30) 

The circle shows the bifurcation value L 2 = 5/2 of the n = m = 1 type-Ill solution 
from the N — 1 type-I solution. The stars indicate the three bifurcating values 
(L 2 = 15/2, 10, 23/2) of the n = 2, m = 1, 2, 3 type-Ill solution from the N = 2 type- 
I solution (energies of the n = 2 type-Ill solutions are not plotted). The numbers in 
parentheses in the figure represent the number of negative modes and the number of 
zero modes respectively of the corresponding solution. They follow from the stability 
analysis which is the content of the next section. 

3 Stability 

3.1 Type-I solutions 

We start with the stability analysis of the type-I solutions @. This set contains the 
lowest energy non-trivial solution, the branch W\ of [JTJ] , which was shown in to be 
for L > \j5/2 classically stable and therefore to be a soliton of the model. Do there 
exist more solitons in this class of solutions? 

To analyse the stability of these solutions we will adopt the point of view of 
hidden algebra and of differential operators preserving finite dimensional spaces of 
polynomials ||, which in the case at hand is a rather straightforward application of 
Fourier analysis. 

Perturbing as usual the fields <p a , a — 1, 2 around the classical solution (§), denoted 
here by 0^, 

(j) a {x) = (j)^(x) + rj a (x) exjp(-iut) (31) 
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leads to the following equation for the normal modes : 



A(N,L) \=uj 2 \ (32) 




N 2 

MN,L) = ~ + 2(1-^ ) " I --j (33) 

sc s 2 ' 

where u 2 is the eigenvalue and, to simplify notation, we posed c = cos(Nx/L) and 
s = sin(Nx/L). 

The complete list of eigenvalues of the operator A(N, L) for N > 1 can be ob- 
tained by classifying its invariant subspaces. This can be done by using the Fourier 
decomposition. After some algebra, one can check that for an integer n > N the 
following finite dimensional vector spaces are preserved by A(N, L) 

px px 

V n = Span{a p cos— + [3 p sm — , p — n = (mod 2N), \p\ < n} , 

Lj Lj 

~ px px 

V n = Span{a p sin— — j3 p cos— , p — n = (mod 2N), \p\ < n} (34) 

Lj Lj 

provided 

a p = P P if p - n + 2N > (35) 

while the other constants a pi (3 p are arbitrary. The operator A(N,L) can then be 
diagonalized on each of the finite dimensional vector space above, leading to a set of 
algebraic equations. One finds finally that the eigenvalues of A(N, L) on V n read 

lu 2 (N, k, ±1) = ((L 2 - N 2 ) +k 2 ± J(L 2 - N 2 ) 2 + 4iV 2 £; 2 ) (36) 

Lj ^ ' 

for k — 0, 1, 2, . . . and similarly on V n for k = 1, 2, 3, . . .. The spectrum given by (|5(f) 
fits exactly with ((^) in the limit L = N for all N i.e. when L approaches the points 
of bifurcation of the type-I solutions from (|]). 

We can now discuss the stability of the solutions of type-I. First remark that 
ui 2 (N, 0, —1) = 0, it corresponds to the zero mode related to the invariance of the equa- 
tions under translations of x. For N fixed and L slightly greater than N, the solution 
(H) possess 4N — 2 negative modes corresponding to e = — 1 and k — 1, 2, 3, ... , 2N — 1 

9 



in (^), remembering that they are twice degenerate. When L increases, more and 
more of these eigenmodes become positive, crossing zero at the values 

L 2 = ^{6N 2 -m 2 ) , m = 1,2,...2JV-1 , (37) 

i.e. (cf. (0)) at those values of L where the solutions of type-Ill with n = N bifurcate 
from the solution of type-I. For 

L 2 >L 2 r (iV)^(6iV 2 -l) (38) 

all the modes are positive and (§) are classically stable solitons. These results are 
illustrated on Fig. [2] for N = 1 and N = 2. The numbers in parenthesis represent 
the number of negative and of zero modes of the corresponding branch. 

3.2 Alternative approach for Type-I 

The stability analysis could also be performed by considering the polar decomposition 
of the Higgs fields @ 

0i = Fcos9 , 2 = Fsin6 (39) 

The quadratic operator Q associated with this paramerization of the fields admits 
the following invariant subspaces 

U n = span{— L=(cos^,0),— L=(0,sin^)} (40) 
VvrL L v7rL Li 



U n = span{ — ^^=(sin — p-, 0), — ^^=(0, cos (41) 



which, by Fourier theory, cover the whole relevant Hilbert space of periodic functions. 
In the normalized basis ([10]), the eigenvalues of Q on the subspace U n read 

u 2 (N,n,±l) = -L{n 2 L 2 + A 2 (n 2 + 2L 2 ) 



± ^/A 4 (n 4 - 20n 2 L 2 + 4L 4 ) + 2n 2 L 2 A 2 (WL 2 - n 2 ) + n 4 L 4 } (42) 

where we defined A 2 = L 2 — N 2 . The eigenvalues corresponding to U n are identical 
to (0). 

10 



These values fail to approach the uj 2 of Eq.(^) in the limit L — > N. This is due 
to the fact that the parametrization (|39D is singular about the solution 0i = 02 = 0. 
Remarkably though, the values (^) and (f42|) have zero crossing at exactly the same 
values of L 2 and the conclusions about the stability of the solutions are identical in 
the two approaches. 

3.3 Type-II and III solutions 

The stability analysis about the solution of type-II can immediately be carried out. 
The equations for the fluctuations about the solution fllpp decouple to take the form 
of the Lame equations: 

{ -£ x + QkSn\y,k)}r ]l = nfo (43) 
{ --^ + 2fc 2 sn 2 (y,fc)}r /2 =nfo (44) 

where we posed 

y = v /TTPx , nl = (ujl + l)(k 2 + l) , a =1,2 (45) 

u) 2 being the effective eigenvalue of the relevant operator. Equations ( fI3j ) and (fP|) 
admit five and three algebraic modes respectively, with corresponding eigenvalues 

Ql : 4 + k 2 , 1 + 4k 2 , 1 + k 2 , 2(1 + k 2 ) ± 2Vl - k 2 + k 4 (46) 

fi 2 : l + k 2 ,l,k 2 (47) 

The corresponding values of uj 2 follow immediately from (fE^); they have signature 
(+, +,0,+,— ) and (0,—,—) respectively for (fjBp and (K7|). Each of the equations 
(El)'©) therefore leads to a zero mode of the solution (|T0|), their origin was dis- 
cussed in Sect. 2. It is a property of the Lame equation that the solutions determined 
algebraically correspond to the solutions of lowest eigenvalues. The remaining part 
of the spectrum therefore consists of positive eigenmodes. The spectrum of Eq. (|43f) 
was studied perturbatively in |4j], while the relation between the Lame equation and 
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the Manton-Samols sphalerons was first pointed out in ||. The presence of nega- 
tive modes in the small oscillation spectrum of all type-II solutions means that no 
classically stable soliton exists among them. 

The stability equation associated with the type-Ill solutions seems to be more 
difficult, mainly due to the lack of analytical expression for Q(x) in Eq. @. We have 
not obtained convincing algebraic expressions for their normal modes but a few of 
their properties can be pointed out. First, these solutions should possess a double 
zero mode due to the invariance of the solution (0), (|l5l) under translations of x and 
internal 0(2) rotation. Likely the eigenvalues relative to these solutions connect to the 
values fl3"o| ) in the limit L 2 — > (6n 2 — m 2 )/2, N = n. We conjecture that the solutions 
of type-Ill always possess at least one negative mode. In the case n — m — 1, the 
negative mode should be unique and meet the two zero modes at the point L = J5/2, 
l) 2 = 0, indicated by the star of Fig. [3]. 

4 Conclusion 

The model ([I]) on S 1 admits a rich set of solutions which bifurcate from each other at 
critical values of the spatial radius L. A detailed analytical study of these solutions 
was presented. Their stability analysis was carried out, all classically stable solitons 
were identified, together with the range of the parameter L for which they are stable. 
The bifurcation pattern contains as a subset the branches corresponding to the solu- 
tions found in Q , and it is found to be richer than that of the gauged Higgs model 
||, the gauged version of (0). 

The spectrum of modes of small fluctuations around a classical solution is shown 
to depend on the parametrization employed for the complex scalar. In particular, we 
explicitly demonstrated this fact for the normal modes around the type-I solutions. 
As one varies the parameter L to approach the bifurcation point L = N along the 
type-I solution with winding N, the correct spectrum is the one which coincides with 
that of the main branch (^). 
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In this work we concentrated on static solutions. The study of time dependent 
ones, to search for the analogs of the breather solutions of the sine-Gordon equation, 
is particularly interesting. They may well be the prototypes of stable "breathing 
membranes" in realistic particle physics models. 
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Figure Captions 



• Figure 1. The solutions of eqs. (|18|) , ([19]) are plotted as functions of k 2 for 
two values of n/m. The dashed line indicates the limit (PI). 

• Figure 2. The energies of some of the solutions are plotted as functions of 
the parameter L. The circle indicates the bifurcation point (L 2 = 5/2) on 
the N = 1 type-I solution. The stars indicate the three bifurcation points 
[1? = 15/2,10,23/2) on the N = 2 type-I solution. The two numbers in 
parentheses refer to the number of negative modes and of zero modes of the 
corresponding solution. 

• Figure 3. The values (|^) are plotted as functions of L for j = 0,1,2 (solid lines), 
together with the values u 2 {l, 1, -1), lu 2 (1, 0, -1), u 2 (l, 0, 1), cu 2 (l, 2, -1), 
uj 2 (1, 1, 1), co> 2 (l,3, —1) of ( p6|) (dotted lines). The numbers indicate the multi- 
plicity of the eigenvalues. 
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